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lized  describing  function  is  derived  for  a  dissymmetric 
relay  with  D.  C,  co.iponencs  in  both  the  input  and  output  paths,  i-lethods 
of  evaluating  the  D.  C.  components  are  discussed.  The  describing  func- 
tion is  applied  to  the  analysis  of  step  and  disturbance  inputs. 

The  describing  function  is  ...odified  to  allow  computation  of  the 
D.  C.  and  sinusoidal  components  of  the  error  signal  when  a  ramp  input 
is  applied.   The  applicability  of  the  method  in  various  situations  is 
discussed. 

The  validity  of  the  describing  function  analysis  xvas  investigated 
by   comparing  calculated  responses  with  actual  system  response,  as  simu- 
lated on  the  analog  computer. 

The  authors  wish  to  express  their  appreciation  for  the  invaluable 
counsel  and  perceptive  suggestions  provided  by  Professor  G.  J.  Thaler  as 
faculty  advisor  for  this  project. 
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.  .  SRI 
:.  PI  LODUCTK 

In  recent  years,  the  describing  function  method  has  been  widely 
used  in  the  analysis  of  servome onanisms  containing  non-linear  elements 
(1).   The  basic  describing  function  technique  consists  of  "linearising" 
the  non-linear  portion  of  the  system  by  assuming  that: 

1 .  The  input  to  the  non-linear  device  is  a  pure  sine  wave  of 

known  amplitude. 

2.  The  higher  harmonics  in  the  output  of  the  non-linear  device 

/  be  neglected. 

.-  representative  describing  function  is  then  calculated  by  com- 
puting the  a  plitude  gain  and  phase  shift  between  the  input  and  the 
fundamental  frequency  of  the  fourier  series  expansion  of  the  output 
wave!   These  techniques  are  thoroughly  discussed  in  standard  texts  (2). 

In  some  systems,  however,  the  presence  of  the  non-linearity 
may  cause  a  D.  Ce  component  to  exist  in  either  the  input  or  output 

bh,  or  both.  This  is  particularly  true  with  dissymmetric  non- 
linearities  (3).  A  similar  situation  can  be  caused  in  many  symmetric 
systems  by  introducing  a  ramp  input.  In  either  case,  the  D.  C„  com- 
ponent may  cause  a  change  in  the  magnitude  and/or  phase  of  the  funda- 
mental frequency  of  the  output  from  the  non-linear  element,  and  it 
may  result  in  a  significant  change  in  the  stability  characteristics 
of  the  syste  i. 

It  is  therefore  advantageous  to  define  a  generalized  describing 
function,  based  on  the  following  assumptions: 


1 ,  linear  element  is  a  pure  sine  xrave 
.  -  .  C.  component. 

2.  The  higher  harmonics  in  the  output  of  the  non-linear 

element  :ay  be  neglected.  In  this  case,  the  output  of 
the  non-linear  element  .nay  be  reduced  to  a  pure  sine 
wave  superimposed  on  a  pure  J.  C.  component,  where 
both  terms  may  be  evaluated  from  the  fourier  series 
expansion  of  the  output  wave. 
The  generalized  describing  function  still  consists  of  the 
amplitude  gain  and  phase  shift  between  the  sinusoidal  component 
of  the  input  and  the  fundamental  frequency  component  of  the  output 
wave . 

The  purpose  of  this  report  is  to  develop  the  generalized  de- 
scribing function  for  a  dissymmetric  relay  and  to  investigate  the 
characteristics  of  servome onanisms  containing  such  elements.  The 
report  includes  a  discussion  of  methods  of  determining  the  magnitudes 
of  the  D,  C,  components;  an  analysis  of  system  stability  for  step  or 
disturbance  inputs;  and  an  analysis  and  discussion  of  the  system 
response  to  step  and  ramp  inputs. 


.  2R  II 
mi  of  d.  c.  sig:;al  levels 

2. 1   General  Considerations 

In  all  cases  investigated  in  this  report,  it  will  be  assumed 
that  the  block  dia grain  of  the  servomechanism  can  be  reduced  to  the 
for:::  shown  in  figure  2.1,  there  G*(s)  and  G?(s)  represent  the  linear 
portions  of  a  type  1  servomechanism,  and  G,  is  the  linearized  describing 
function  -which  represents  the  non-linear  element  in  the  system. 


a 

Kt) 

o(t) 

©i 

G^s) 

ad 

a2(s) 

©0 

Figure  2.1  Block  diagram  of  a  servomechanism  containing  a  dissymmetric 
non-linearity. 


In  addition,  it  will  be  assumed  that  the  servo  loop  acts  suf- 
ficiently like  a  low  pass  filter  so  that  only   the  first  harmonic  and 
the  D.  G.  component  of  0(t)  are  transmitted  around  the  loop.  There- 
fore, for  steady  state  operation  with  sinusoidal,  step,  ramp,  or  dis- 
turbance inputs,  we  can  approximate  the  input  to  the  non-linearity 
with  equation  2, 1 ;  and  the  effective  portion  of  the  output  from  the 
non-linearity  can  be  approximated  by  either  equation  2,2  or  by  squa- 
tions  2,3  through  2.5. 


2.1       I(t)  =  vi  +  Xsin  ut 


2.2  0(t)  =  ~°-  +  a-cosut  +  b^sinwt 

2  '  n 

2.3  O(t)  =  ^°-  +  Csin(u>t  +  ^) 


2.4  c  =  yai2  +  b.,2 


2.5      /  =  tan-1       21 


It  will  also  be  assumed  that  G-(s)  represents  a  compensator 
having  a  transfer  function  of  the  for:  expressed  in  equation  2,6;  and 
that  Qp(s)  represents  a  r.otor-load  combination  having  a  transfer  func- 
tion of  the  fori  expressed  in  equation  2.7. 

2.6  Gl(s)  =k,  liHiil 

TT  (s  +  ?i) 

where     z.   ^  0  and  p^_  ^  0 

2.7  0,(.)  =  k,    "<■'  +  zi\ 

where     z .-,  j£  0  and  p  .  /  0 

The  basic  Laplace  transform  equations  for  such  a  system  are: 
2.3       £    =    ©i  -    ©0 


2.9       Kt)   =  G^s)^  -©0) 


2.10      eo  =  G-2(s)0(t) 


0   2 


Df  C,  Levels  in  the  Cutoab  Path 

In  the  case  of  s':   ,  listurbance,  ramp,  or  sinusoidal  inputs, 
the  D.  C,  component  of  the  output  fro:  the  non-linear  element  depends 
only  on  the  linear  sain  of  (^(s)  and  the  applied  input  function,  and 
is  relatively  easy  to  evaluate. 

The  equation  for  the  Laplace  transform  0(t)  can  be  obtained 
by  using  only  the  transfer  functions  of  the  linear  portion  of  the 
system  as  in  equation  2,11.  Since  the  steady  state  sinusoidal  com- 
ponents have  no  effect  on  the  D,  C.  components  in  the  linear  portion  of 
the  system,  the  principle  of  superposition  can  be  used  to  eliminate  the 
transforms  of  these  components  when  solving  for  the  D.  C.  components 
only.  The  final  value  theorem  can  then  be  applied  to  obtain  equation 


2.12,  where  C(t)^,  ©j^,  and  I(t)^  represent  the  Laplace  transforms  of 
the  non-sinusoidal  components  of  0(t),  ^,  and  I(t),  respectively. 

<3o  _  1 


2.11   0(t)  = 


TS7  "  uoXsj 


1     w 


2.12  t°   =  sO(t) 

9 


d 


cQ^T 


s=o 


6  J^d 

id   ^JTsT 


s=o 


In  order  to  simplify  the  following  calculations,  the  limits  of 
G.j(s)  and  G?(s)  as  s  approaches  zero  can  be  partially  evaluated  to 
get  equations  2.13  a^d  2.1^.  These  results  can  then  be  sul   '  iterl 
into  equation  2.12  to  get  equation  2.15. 

TT(s  +  zj 


2.13 
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s^o 
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s-*o 
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For  a  step  input,  or  disturbance,  of  magnitude  A,  the  Laplace 
transform  of  ^  is  given  by   equation  2.16,  For  a  steady  state  analysis, 
the  initial  time  is  immaterial,  and  consequently  the  Laplace  trans f or n 
of  I{t)j  can  be  represented  by  equation  2.1?  and  the  D.  G.  component  of 
C(t)  can  therefore  be  computed  by  equation  2.1  . 


2.16  9.  =  Au(t)  =  £ 

1  s 

2.17  Kt)d=| 


2.1   Sg 
2 


''Y2 


v1 


=  0 


J  s^.o 


A  sinilar  argument  can  be  applied  to  sinusoidal  inputs,  since 


in  this  case 


,  &, 


id 


xs  zero. 


For  ramp  inputs  of  magnitude  3,  the  Laplace  transform  of  the 
input  is  given  by  equation  2.19,  and  the  D.  C.  component  of  0(t)  can 
be  evaluated  by  equation  2.20. 


2.19  e}  =  Btu(t)  = 


2.20  ^-  =  r£ 


-2 1 


^2   s  k 


V1 


s-^o 


For  an  ncn  order  s        le  s  term  in  the  denominator  of 
G£(s)  would  be  changed  to  sn,   Equation  2*21  can  be  derived  for  such 
a  system  by  following  a  line  of  reasoning  similar  to  that  used  for  the 
first  order  system.  Evaluation  of  equation  2„21  shovre  that  the  d.c, 
component  of  0(t)  will  be  zero  for  both  step  and  ramp  inputs  for  any 
system  of  order  higher  than  ones 


o  01   aO  _  Sn 


wid  -  sir- 
V1 


s*.  0 


2.3  D9  C.  Levels  in  the  Inout  Path 

In  the  preceding  discussion,  there  was  no  limitation  on  the 
nature  of  the  non-linearity  represented  by  G^,  since  the  D.  C,  com- 
ponent in  the  output  path  is  independent  of  G^»  In  general,  however, 
the  D,  C,  level  in  the  input  path  will  not  only  depend  on  the  nature 
of  the  system  input  and  the  characteristics  of  the  linear  portion  of 
the  system,  but  will  also  depend  on  the  nature  of  the  non-linear  ele- 
ment . 

In  many  cases  the  solution  for  the  Ds  C,  component  in  the  input 
path  can  be  extremely  complicated,  and  may  result  in  transcendental 
equations  which  can  best  be  solved  by  graphical  means  (3). 

It  should  be  noted,  however,  that  in  many  cases  it  is  unnecessary 
to  solve  explicitly  for  the  Ds  G„  level  in  the  input  path  in  order  to 
obtain  the  describing  function  curve,  This  fact  is  particularly  well 
illustrated  in  the  analysis  of  the  relay  servo  with  a  ramp  input,  given 
in  section  ^.2,  In  other  cases,  such  as  the  analyses  for  step  inputs 
given  sections  3.2  and  3.3»  the  solution  for  the  D.  G.  level  indicates 
that  for  frequency  responses  analysis,  the  D.  C,  level  in  the  input  path 


ls  treated  as  a  sym- 
metrical compone 


ran  in 

;i ..  .  savo 

3.1  Discussion  of  the  Describing  Function. 

..  pendix  I  contains  the  derivation  of  the  generalized  describing 
function  for  a  dissymmetric  relay  with  the  following  characteristics: 

1 .  The  pull-in  voltage  for  a  positive  input  is  not  equal  to 

the  pull-in  voltage  for  a  negative  input. 

2.  The  drop-out  voltage  for  a  positive  input  is  not  equal  to 

the  drop-out  voltage  for  a  negative  input. 

3.  The  hysteresis  for  a  positive  input  is  not  necessarily  equal 

to  the  hysteresis  for  a  negative  input. 
U-,     The  positive  relay  output  is  not  equal  to  the  negative  re- 
lay output. 
The  input  vs.  output  characteristics  of  such  a  relay  are  shown 
in  figure  1-1  of  Appendix  I. 

Since  most  practical  relays  have  positive  and  negative  outputs 
of  equal  magnitudes,  no  attempt  was  made  to  analyze  a  relay  with  un- 
equal outputs.  Therefore,  in  the  remainder  of  this  report,  it  will  be 
assumed  that  the  relay  output  voltages  are  of  equal  magnitude. 

3.2  The  Describing  Function  for  a  Relay  with  no  Hysteresis. 

The  input  vs.  output  characteristics  for  a  dissymmetric  relay 
with  no  hysteresis  can  be  represented  by  figure  3»1. 

Equation  2.13  shows  that  for  a  servome onanism  of  the  type  shown 
in  figure  2.1,  the  D.  C.  component  of  the  relay  output  is  zero  for 
steady  state  operation  after  the  application  of  a  sinusoidal,  step,  or 
disturbance  input.  Combining  this  result  with  equation  1-15  of  Appendix 
I  gives: 
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r 
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u» 


*-A 


u>t  3  7r-^( 


out  *  TT+  3. 
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Figure  3.1  Transfer  characteristics  of  a  dissymmetric  relay  with  no 
hysteresis. 
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.1   ia  =  S(^  -fiz) 

Tf      l 

3.2  B,  =0? 

The  angle  relationships  of  equations  1-9  and  1-11  of  Appendix  I 
can  nou  be  introduced  to  solve  for  the  D.  C.  component  of  the  input 
uave  and  to  solve  for  A  and  ^   as  shown  below. 

d-  -  W        n  do  +  W 


.3  sin«j  =  — L_ =  sinPp  =  £x 


^        ■  2 


3.5  sin^j  =  — L — ~Jk    -   sin ^ 


For  a  relay  with  no  hysteresis  the  expression  for  a..,  as  derived 
in  Appendix  I,  is 

3.6  a.j  =  0 

ie  expression  for  b-  is  given  by  equation  1-16  of  Appendix  I, 
combining  this  equation  with  equation  3.2  gives 

3.7  b*   =  £i(cos&>  +  cos  $,)  =  ilcosA 

7T  77  ' 


■jjp 


/1  -  sin2&    =  s 


d1  •  d2 


The  equation  for  the  generalized  describing  function  is; 


3.8  Ife-i/^+V     /*>_£• 


For  the  dissymmetric  relay  with  no  hysteresis,  the  describing 
function  is: 
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.    G .  =  £ 


d1  +  d2 


for  X  >  d1  *  d2 


Table  I  is  a  list  of  the  values  of  the  generalized  describing 
function  for  a  relay  with  d^  +  d2  =  5.2,  and  1=1.  A  graph  of  this 
describing  function  is  shown  in  figure  3.2. 


X 

Gd 

Gd 

Gd"1    (db.) 

2.6 

0 

CO 

o° 

3.0 

.212 

4.73 

13.5 

3.69 

.257 

3. 

11.5 

4.0 

.242 

4.13 

12.25 

5.0 

.217 

Mo 

13.2 

6.0 

.191 

5.23 

14.3 

7.0 

.169 

5.92 

15.4 

8.0 

.150 

6.65 

16.4 

10.0 

.123 

S.13 

17.2 

... 

Table  I 


Values  of  G^  for  a  relay  with  no  hysteresis. 
d1  +  d2  =  5.2 


Equation  3.9  applies  only  for  cases  in  uhich  the  relay  operates 
on  both  sides  during  each  cycle.  It  is  not  necessary  to  consider  one- 
sided operation  at  this  point,  however,  since  this  type  of  operation 
requires  that  the  D.  C.  level  in  the  relay  output  roust  be  different  from 
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zero.  It  has  previously  been  s]     hat  for  a  servomechanis  i  of  the 
type  described  in  Chapter  II,  the  stead/  state  D,  C,  component  in  the 
relay  output  must  be  zero  for  sinusoidal,  step,  or  disturbance  inputs. 

If  the  D.  C.  level  in  the  output  is  allowed  to  be  different  from 
zero,  as  in  the  case  of  a  ramp  input,  then  both  one-sided  and  two- 
sided  operation  is  possible.  The  subject  of  ra     >uts  is  treated  in 
Chapter  I'.ra 
3.3   'he  DescriM.;  b  .c^ion  for  a  Relay  with  Hysteresis. 

-be  input  vs.  output  characteristics  for  a  dissymmetric  relay 
with  hysteresis  can  be  represented  by  figure  3.3  «  In  this  case,  the 
co  ibined  results  of  equation  2.13  from  Section  2,2,  and  equation  1-12 
fro  i  Appendix  I  gives: 

3.10  f2  =  L-(/9  +<*  -&     -*  )  =  0 

2   2fT  2    2    1    1 

3.11  fiy  ,^=6z    +  * 

Combining  the  angle  relationships  given  by  equations  1-3  through 
1-11  with  equation  1-13  of  Appendix  I,  we  obtain  equation  3.12. 


3.12  &1  =- 


[1  -  Pj  +ib-^2 
X        X 


If  the  hysteresis  is  the  same  for  both  sides  of  the  relay, 
d..  -  p.  =  d2  -  p2,  and  the  angle  relationships  can  be  evaluated  as 
follows : 

3.13  sinA  -  sin  <tf  =  d1  -  ?,1  =  d2  "  p2 


sin^>  -  sino( 


._ 


2 
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ut  »  «(, 


•      tO  t  *  IT  -  &, 


X  a  n  +  °£ 


Ujt*2lT'fi3 


Figure  3.3     Transfer  characteristics  of  a  dissymmetric  relay  with 
hysteresis. 


3.14     sin/S,   -  siiW,   =  2sinK^  -  tf^cosiC^  +  <, ) 

=  2sini(&,  -  o(2)cosi(^2  +^2) 

By  combinin  ;  t  3  results  of  equations  3.11 t  3.13i  and  3.14,  it 
can  be  shown  thatoCi  =  o<2  and  #..  =  /99,     The  D.  C.  component  of  the 

input  can  then  be  evaluated  by  equation  3.16»  and  the  angles  can  be 
evaluated  by  equations  3.17  and  3.13. 


J.  15  sin^l,  =  — i— ■ —  =  -=-» —  =  sin  4, 


3.16  -  =  a1  "  Cl2  _  ?1  -  ?2 

3.17  sin*.  =  sjjioc  =  P1  ""  P2  =  R 

1      *•     ?.x 


o 


US  sin^  =sin/|=iLH 


2      2X      d 

Ehe  general  expression  for  b..  is  given  by   equation  1-14  of 
Appendix  I.  Evaluating  this  equation  by  inserting  the  results  from 
equations  3.17  and  3.18  gives  equation  3.19.  The  describing  function 
equation,  3.20,  is  then  obtained  by  substituting  equations  3.12  and 
3.19  into  equation  3.2.  The  describing  function  in  this  instance  also 
reduces  to  exactly  the  sane  form  as  the  describing  function  for  the  cor- 
responding  symmetric  relay. 

3.19  b.  =  —  (cos£  +  cosO 
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3.20  ^y~~/  +  (/77~^/7T?  y 


i>  =   tan-1        Tid   "  *p 


Hd2  +  1  -  R32 


for  X  >  Pj  +  Pg 
2 


Table  II  contains  the  values  of  the  generalized  describing  func- 
tion for  a  relay  with  d..  +  d0  =  6   and  p1  +  pg  =  10,  A  graph  of  this 
describing  function  is  shown  in  figure  3.';-. 
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X 

■d 

v1 

V1(db) 

Phase  Angle 

; 

.1135 

-. 

18. 

26.6 

5.2 

.1^13 

7.08 

17.0 

19.5 

5.5 

.1505 

6t65 

16.4 

16.2 

6 

.15^3 

6.46 

16.2 

13.2 

7 

.1  I 

o.75 

16.7 

10.1 

0 

.1372 

7.23 

17.2 

8.33 

9 

.1265 

7.9 

17.9 

7.15 

10 

.11 

3.56 

13.7 

0.26 

11 

.1 

9.27 

19.3 

5.62 

12 

.100 

10.0 

20 

5.1 

15 

.03018 

12.22 

21.7 

3.97 

20 

.0624 

16.0 

24. 1 

2.92 

30 

.0422 

23.7 

27.4 

1.92 

40 

.0316 

31.6 

29.9 

1.44 

50 

.0254 

39.35 

31.8 

1.15 

Table  II    Values  of  Gj  for  a  relay  with  hysteresis. 
d1  +  d2  =  6  and  P1  +  P2  =  10 
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3  same  for  both  sides  of  the  relay, 
the  transcendental  equations  which  must  be  solved  in  order  to  determine 
the  D,  C,  component  of  I(t)  become  much  more  complex.  They  are  difficult 
to  solve  with  trigonometric  relationships  alone,  but  it  is  relatively 
easy  to  solve  them  by  graphical  methods  using  the  construction  shown  in 
figure  3.'  >       'his  construction  can  be  developed  in  the  following  manner: 

1 .  Selected  a  protractor  with  radius  R. 

2.  Calculate  the  following  quantities  for  each  value  of  X  to  be 

investigated, 

X 

x  x 

3.  Construct  figure  3.5  using  the  calculated  dimensions. 
k-t     Move  the  protractor  along  the  vertical  axis  until 

1  +  °S  =  2  +  °^2  as  s^01Jn  on  fiSure  3.6.  This  satisfies 
equation  3.11. 
The  construction  guarantees  that  equations  1-8  through  1-11  are 
satisfied  if  and  only  **«. ,  ft  ,   °C  and  P~   are  selected  as  shown  in  figure 
3.6,  and  the  D,  C,  component  of  I(t)  is  established  by  the  center  of  the 
protractor.  Figure  3.6  therefore  represents  the  desired  solution  for  the 
angles,  and  ;.  =  §L 

The  generalized  describing  function  can  then  be  evaluated  using 
equations  3.S.  3.12,  and  1-1^.  Table  III  gives  the  values  for  such  a 
describing  function  for  a  relay  with  p1  -  3»  p?  =  7»  cL  =  2  and  d^  =  LK 
s  describing  function  is  plotted  in  figure  3.7. 
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V 

- 1 

, 

u  ? 

0 

i * 

Figure.  3.5 


Figure  3.6 


Figures  3.5  and  3,6     G^arhical  construction  for  solution  of  equations 

3.11  and  1-3  through' 1-11. 
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- 

■ 

4 

Gd"1  (do) 

. 

-1.800 

.1290 

21.4 

17.3 

6.0 

-1.615 

.1575 

13.45 

16.25 

7.0 

-1.562 

.1525 

10.2 

16.25 

8.0 

-1.540 

.1370 

S.35 

17.1 

10.0 

-1.520 

.1164 

6.3 

10. 7 

12.0 

! 

-1.520 

.0999 

5.1 

20.0 

15.0 

-1.510 

.on 

4.0 

21.7 

Mote: 

Values  for  G^ 

and 

t 

for  la:- 

value 

s  of  X 

are 

approximately 

• 

the 

sane  as  corresponding  values 

in 

table 

II 

Table  III    Values  of  the  generalized  rela;/  describing  function  for 

p1  =  ^»  p2  =  7»  d1  =  2»  and  d2  =  lv* 
3.4  Cxroarison  of  Sy  netric  and  Diss  -metric  Relays. 

Equations  3.8  and  3.20  indicate  that  the  generalized  describing 
functions  for  dissy  metric  relays  are  exactly  the  sane  as  the  corresponding 
describing  functions  for  symmetric  relays,  provided  the  dissymmetric  relays 
have  equal  hysteresis  on  each  side.  A  servomechardsm  containing  such  a 
dissymmetric  relay  would  therefore  respond  to  a  sinusoidal,  step,  or 
disturbance  input  in  a  manner  quite  similar  to  a  servo  containing  a  sym- 
metric relay. 

If  the  servo  response  results  in  a  limit  cycle  or  sinusoidal  out- 
put, the  output  oscillations  for  the  symmetrical  relay  would  be  centered 
about  the  position  commanded  by  the  D.  C.  portion  of  ©. ,  while  the  output 
oscillations  for  the  dissymmetrical  relay  would  be  offset  slightly  to 
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create  the  necessary  J,  C.  Tec-.  re1:  signal  at  the  input  of  the  relay. 

Thus  the  effect  of  die      ry  is  to  introduce  a  steady  state  D,  C.  error, 

If  the  dissymmetric  relay  does  not  have  equal  hysteresis  zones  on 
each  side,  the  generalized  describing  function  a  -  be  quite  different 
fro:,:  the  describing  function  which  uould  be  obtained  by  assuming  equal 
hysteresis  zones  and  using  equation  3.20,  A  comparison  of  figures  3.^ 
and  3.7  indicates  that  these  differences  will  be  relatively  minor  in  most 
practical  cases,  and  that  adequate  results  can  be  obtained  by  using  equa- 
tion 3,20, 
3.5  Verification  of  Proposed  Describing  Functions, 

In  order  to  verify  the  describing  functions  previously  derived,  a 
position  control  system  consisting  of  the  dissymmetric  relay  and  third 
order  linear  components  was  postulated.  Linear  transfer  functions  with 
typically  low  pass  characteristics  were  selected  for  this  system. 

Figure  3*3  is  a  schematic  of  the  control  system  together  with  the 
analog  computer  simulation.  This  simulation  provides  for  operation  of 
the  relay  under  varying  degrees  of  dissymmetry  and  hysteresis.   (-'.'-) 
Figure  3.9  is  a  Nichols  Plot  of  the  linear  transfer  functions  with  an 
arbitrary  gain  factor  of  20, 

Case  I  Dissymmetric  Relay  '.."ithout  Hysteresis. 

The  relay  was  operated  without  hysteresis  and  with  a  constant 
total  dead  zone  of  5.1  volts.  The  plus  and  minus  pull-in  voltages  were 
varied  on  a  series  of  runs  to  obtain  various  conditions  of  dissymmetry. 
Step  inputs  of  10  volts  and  15  volts  were  applied  to  the  system.  The 
gain  for  the  linear  portion  of  the  system  was  adjusted  to  different 
values  so  as  to  place  the  linear  curve  in  the  desired  position  on  the 
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Nichols  Plot,  thus  pro        series  of  limit  cycles  of  various 
amplitudes. 

fable  IV  is  a  summary  of  the  results  obtained.  Typical  Brush  re- 
cordings for  several  values  of  system  gain  and  input  magnitudes  are  dis- 
played in  figures  3.10  through  3*12,  which  illustrate  the  predicted  steady 
state  d.c,  error  levels  and  sinusoidal  characteristics. 


System  Observed   Limit  Cycle  Predicted  Limit  Cycle 


Run 

Relay 
Pull-in 

3yst 
Gai: 

1 

+2.o,-2.5 

88 

2 

+3.1,-2.0 

88 

3 

+2.4,-2.7 

132 

4 

44.0,-1.1 

132 

5 

+2.0,-3.1 

132 

6 

+3.5.-1.6 

132 

7 

+2.7,-2.4 

81 

8 

+2.7,-2.4 

63 

>litude 

Frequency 

•  -  iplitude 

Frequency 

4.6 

2 

4.7 

2 

4.6 

2 

4.7 

2 

8,0 

2 

0.0 

2 

7.9 

2 

G.O 

2 

8.1 

2 

8.0 

2 

3.0 

2 

8.0 

2 

4.0 

2 

3.8 

2 

None 

mm 

Rone 

^ 

Table  ]  ,   Experimental  Results  for  Dissymmetric  Relay  Without  Hysteresis. 

Case  II  Dissymmetric  Relay  with  Hysteresis, 

2  relay  was  operated  with  the  following  characteristics:  pull 

in,  '■   .  p  volts  and  -   .  >  volts;  drop  out,  +  4.0  volts  and  -  1.95  volts. 

eral  procedure  was  followed  as  in  Case  I.  <   results 

is  shown  in  table  ,    pical  Brush  recordi  s  -.re  displayed  in  figures 

3.13  thro-    .  i  . 

..  comparison  between  the  predicted  and  observed  values  shown  in 

table  V  indicates  the  validity  of  the  proposed  describing  function  for  a 
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relay  with  hysteresis,  11  '.  that  a  snail  difference  between 

the  observed  :     'edicted  values  is  to  be  expected  since  the  experimental 
relay  characteristics  were  not  exactly  matched  to  those  used  for  predic- 
tions.  Ihe  predicted  values  were  based  on  figures  J,k  where  p^  +  ?2  =  10 
and  d^  +  Co  -   6.   Ihe  characteristics  of  the  experimental  relay  xrere: 
P1  +  p2  =  9.3  and  d..  +  d0  =  5. .  . 


Run 

Gain 

Observed 
A  rcplitude 

Li  lit  Cycle 
Frequency 

Predicted 
Amplitude 

Limit  Cycle 
Frequency 

1 

620 

39.6 

1. 

kO 

1.93 

2 

29.0 

1.88 

30 

1.90 

3 

300 

19.35 

1.75 

20 

1.35 

4 

172 

12,75 

1.72 

12 

1.75 

5 

7.15 

1.50 

6 

1.51 

rable  V.   [experimental  Results  for  Dissymmetric  Relay  '..ith  Hysteresis. 
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IV 
;   .    1  FUTICTIOrJ  ?OR  A  RAIT  INPUT 

■ • 1  General  Considerations. 

sn  a  relay  servo  is  subjected  to  a  step  or  disturbance  input, 
the  servo  T.dll  either  co.ie  to  rest  uhile  the  error  remains  within  the 
limits  of  the  dead  zone,  or  it  will  enter  a  limit  cycle  which  operates 
each  side  of  the  relay  during  each  cycle.  Similarly,  when  a  sine  input 
is  applied,  the  steady  state  error  signal  will  either  vary  within  the 
limits  of  the  dead  zone,  or  it  will  cause  the  relay  to  operate  on  both 
sides  during  each  cycle, 

.  :-n  subjected  to  a  ramp  input,  however,  the  error  will  never  re- 
-  ain  in  the  dead  zone.  The  servo  output  nay  cone  to  rest,  but  the  in- 
creasi-.  ;  ra       •.!  in  the  input  channel  will  always  cause  the  relay 
to  actuate  again.   Therefore,  the  servo  must  enter  a  limit  cycle  when- 
ever a  ramp  input  is  applied,  unless  the  ran  signal  is  so  large  that 
the  relay  does  not  drop  out.  This  limit  cycle  may  involve  actuation 
of  both  sides  of  the  relay,  or  it  may  involve  or.?      repeated  operation 
of  one  side. 

-■   case  in  which  the  relay  operates  on  only  one  side  is  the 
si  plest,  and  will  be  discussed  in  Section  ^.2.   ]  on  the  relay  operates 
on  both  sides,  the  solution  must  be  obtained  by  Graphical  methods,  as 
explained  in  Section  ^-.3. 
h%2.     Describi:'  ;  .  .  .iction  for  One-sided  Operation. 

he  relay  operates  on  only  one  side,  the  input  vs.  output 
characteristics  can  be  portrayed  '   fi.  ;ure  'J,   .  .e  fourier  expansion  of 
the  output  T.:ave,  which  is  derived  in  Appendix  I,  becomes 
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Figure  b;i     Transfer  characteristics  of  relay  operating  on  only  one  side, 
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'+.1     a0  =  I  IV-   (oC-!-/3)l 
7T  L  J 

i   =  L  (siri/fi  -  sino( ) 
TT 

';._   .j  =  I  (cos£  +  cos  oO 
ten  the  servo  has  s.  ramp  input,  O.  =   'Jtu(t),  the  D.  C.  component 
of  the  rela:r  output  wave  is  given  by  equation  2.20.   rhis  equation  can 
be  combined  with  equation  fy.1  to  get  equation  ^.4 


2.20  fft  =  JL 

2   ^~ 


The  equations  for  the  generalized  describing  function  can  now  be 
obtained  in  the  following  manner,  without  solving  explicitly  for  the  D. 
G.  component  of  the  relay  input. 


^•5    a-j     +  ' 


1 


■:- 


'(sin2/9  +  sin  o<  -  2sin$  sinol) 
(cos  >6  +  cos  <<  +  2cos/}cosc() 


"   -l... 


•77 


(2  +  2coso(cos/Q  -  2sino(sin/8) 


=  412 


=  2 


1  +  cos(<\  +  £  ) 
[If     [1   -  cos^-*<  "*)" 

fft 


1   -  cos 


23  TT 

"to 


JJ..6     iL  =  s^-  sinoC   =  tan  1(^.^0 
b1       cos£  +  coS(?(  2 

*.7     4  =  ban  -1  tl  =  k*-<) 
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At  this  point,  it  is  convenient  to  introduce  the  trigonometric 
identity  given  by  equation  4. 

:,    sin/8-  sino(=  2sin  l(^-o()cos  l/£  +  <^\ 

2         ^ 


=  2  sin  ^  cos 


TT 


B7T 


2  "  Yfr 


.        77 
=  2  sin  ©  sin  -zJL. 

YKV 
v2 


4.9   sin  I   =  sin<9-  sjji<* 
2sin  3^: 
lhV2 

4-.  10  sin  i   = 


2Xsin  -ML 
TK__ 
v2 


The  generalized  describing  function  for  the  ramp  input  during 
one-sided  operation  becomes 


The  relay  will  operate  on  only  one  side  as  long  as  X  is  less 
than  a  critical  value,  X  ,  defined  by  equation  4.12. 

4.12  X  =  SLlS 

c    1  +  sin  oC 

If  the  relay  has  no  hysteresis,  there  should  be  no  discontinuity 
in  the  describing  function  as  the  servo  shifts  from  one-sided  operation 
to  two-sided  operation.  On  the  other  hand,  if  the  relay  has  a  signifi- 
cant amount  of  hysteresis,  there  will,  in  general,  be  a  discontinuity 

when  the  shift  occurs. 
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The  expression  f  .  magnitude  of  Gj  contains  no  minimum  limit 
for  the  amplitude  of  the  sinusoidal  input*  The  phase  angle  expression, 
however,  indicates  that  the  describing  function  is  valid  only  if 

2sin-BE 

Mote  that  if  the  denominator  of  this  expression  is  equal  to  its  inaximu  i 
value,  2,  this  restriction  merely  requires  that  the  peak  to  peak  oscil- 
lations of  the  sine  wave  must  extend  from  one  side  of  the  hysteresis  zone 
to  the  other. 

For  relays  with  absolutely  no  hysteresis,  the  describing  function 
predicts  absolutely  no  phase  shift  due  to  the  relay.  However,  if  the 
relay  contains  at  least  an  infinitesimal  amount  of  hysteresis,  as  all 
real  relays  do,  there  will  be  some  small  value  of  X  for  which  the  magni- 
tude of  Gj  is  defined,  and  for  which  the  phase  shift  will  be  90  degrees. 
Consequently,  the  effect  of  hysteresis  must  be  considered  to  be  much 
.■lore  significant  in  the  case  of  ramp  inputs  than  for  other  types  of  inputs. 

The  describing  function  defined  by  equation  4.11  can  be  plotted 
in  non-dimensionalized  form  by  defining  the  following  quantities 

*.»  g^  2i<£LLi> 


1    Z(\   -  cos2-E. 
^  J       \  2KV2 


4.15   sin  0  =  =1 


2Zsin  —M 
*2 
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\.1    3(s)  =  1 :^(s)G2(s) 


?1  -  dt 


4.1?  Dh  =  f  „i  "'"  '??■ 

-  d1 


The  critical  equation  for  stability  beco  :es 

4.13  3dCri(s)32(s)  :  CTdn-(s)  =  -1 

1 


The  plot  of 


n 


,.  /i       can  then  be  used  as  the  ramp  describing 
do* — 


function  for  one-sided  operation  of  any  relay  by  simply  adjusting  the 

gain  of  the  linear  plot  by  the  factor  „J .  A  graph  of  the  de- 

P1  -  d, 

scribing  function  of  equation  4.14  is  plotted  in  figure  4.2. 

maximum  value  of  Z  for  which  the  describing  .function  will  be 
valid  is 

4  19   Z  =    Xe   .  =  ,     p1  +'?Z 

c   p1  -  d.,   Tp^pTcl'^d  -TsZrnjJ 

nc 


1  +  sin  p( 

e  describing  function  plot  can  also  contain  loci  of  equal 
values  of  D^  =  Z(l  +  sinoO.  These  loci  will  pass  through  points  which 
correspond  to  the  ma::;     ralue  of  Z  for  which  one-sided  operation  will 
occur  for  a  relay  with  the  corresponding  ratio  of  dead  zone  to  hysteresis. 
In  order  to  obtain  these  loci,  it  is  most  convenient  to  first  plot 

Z  vs  Z(1  +  sin  o()  for  various  ratios  of  , .,  •?  ,  This  requires  solving 

^r2 

for  the  angle  0  using  equation  4.15  then  solving  for  the  angle  o(   using 
equations  4.4  and  4.7.  The  values  of  Z  can  then  be  obtained  for  constant 
values  of  DTl  and  these  ooints  can  then  be  plotted  on  the  describing  function 
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i 


Figure  4.5     Illustration  of  wave  shape  of  I(t)  in1  non-sinusoidal  region 
of  ranp  describing  function  plot. 
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curve,  lots  of  Z  vs  Z(1  +  s±lio()   for  various 

values  of   3  .   rhe  information  on.  these  curves  was  used  to  obtain  the 

curve  for  D^  =  10,  which  is  shown  on  figure  4.2. 

It  can  also  be  shovni  that  there  is  a  limiting  value  of  phase 

shift,  beyond  which  the  relay  input  can  no  longer  be  sinusoidal.  If 

>  ,5YX  »  sinusoidal  operation  is  possible  only  if 
2 


4.20  -0  S  77-  Jl 


377 


■Uv2 


if  .  ^  ,5^K  ,  sinusoidal  operation  is  possible  only  if 


v2 


4.21  4    <  -JL 

"  Tz 


These  limiting  values  are  also  plotted  on  figure  4,2.  The  mathematical 
proof  of  the  limitation  is  contained  in  Appendix  II, 

The  non-sinusoidal  characteristics  become  more  and  more  pronounced 
as  » —  approaches  zero  or  unity,  and  as  0  approaches  90  degrees.  For 


v2 


example,  for  B  =  .1  Ykv  and  j£  =  81°,  the  solution  of  equations  4.4  and 

2 

4.7  requires  that  o(.  =  172°  and  &  -   0°.  This  requires  a  wave  shape  of 

the  form  illustrated  in  figure  4.5.  The  corresponding  sinusoidal  wave  is 

shown  in  dotted  lines  for  comparison. 

4.3  Experimental  Procedures 

As  a  first  step  toward  the  determination  of  a  describing  function 

for  ramp  inputs,  a  series  of  investigative  runs  were  made  on  analog 

control  systems.  For  this  purpose  two  control  systems  were  simulated  with 

transfer  functions     ^  . ,  and        K         resoectively,  Nichols 

3(3  +  4)      J(j  +  2}(^  +  4) 

Plots  for  these  systems  are  shown  in  figure  4.6. 
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Analysis  of  the  results  of  these  preliminary  runs  together  with 
the  theoretical  concepts  which  had  been  developed  up  to  this  point, 
indicated  the  significance  of  a  parameter  formed  by  the  ratic  .....P...  . 

An  experiment  was  then  designed  to  more  completely  explore  this  feature, 
A  series  of  values  for  the  ratio  > 3  .  ,  ranging  from  0.01  to  0.90, 

YKy2 

were  selected.  Several  runs  were  made  at  each  of  the  selected  values 
of  the  ratio,  by  varying  B  and  K   to  keep  the  ratio  constant.  The 
relay  was  operated  at  a  constant  value  of  hysteresis  and  data  was 
recorded  for  both  the  second  and  third  order  systems.  The  reduced  data 
is  shown  in  table  VI,  and  on  a  Nichols  Plot  in  figure  4.7. 

An  examination  of  the  data  plotted  in  figure  4.7  leads  to  the 
formulation  of  certain  conclusions.  Since  each  set  of  limit  cycle  condi- 
tions, slotted  for  a  constant  value  of    ^  ,  seems  to  lie  on  a  single 

XKV2 

curve  it  appears  that  a  describing  function  curve  exists  for  each  value 
of  this  ratio.  The  fact  that  the  limit  cycle  conditions  obtained  for 
the  second  order  system  substantiate  the  pattern  established  by  the 
third  order  system  even  more  clearly  illustrates  the  existence  of  a 
describing  function.  The  nature  of  these  describing  function  curves  in- 
dicates that  they  will  extend  across  all  phase  angles  from  -90°  to  -130° 
if  the  curves  are  extrapolated  beyond  the  region  of  plotted  data.  This 
agrees  with  intuitive  reasoning  that  a  limit  cycle  must  exist  for  all 
type  I  systems  with  ramp  inputs,  unless  the  system  is  incapable  of  fol- 
lowing the  ramp. 
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YIC 

lily 

- 

Order 

Svste 

Gd"*1 
db 

Second  Order  Svste 
Freq.           X         Phase 
rad/sec     volts     Angle 

3^-1 

ud 

db 

3 

'I. 

rad/sec 

X 

volts 

Phase 
Angle 

.5 

.01 

.924 

2.0 

103 

34.5 

.5 

.02 

.342 

4.45 

109 

36.2 

1.03 

1.47 

106 

26.7 

1.0 

.02 

1.515 

1.94 

111 

29.9 

.5 

.05 

.533 

2.40 

125 

23.5 

1.33 

0.95 

139 

16.5 

1.0 

.05 

.61 

3.67 

130 

23.2 

1.965 

1.21 

116 

19.0 

1.5 

.05 

.603 

4.8 

134 

30.5 

2.40 

1.4 

121 

20.2 

2.0 

.05 

2.73 

1.56 

124 

21.3 

2.5 

.05 

3.14 

1.6? 

128 

21.8 

.5 

.  1 

.6.':- 

1.65 

131 

15.7 

1.65 

.72 

112 

3.9 

1.0 

•  i 

.796 

2.26 

1  ';0 

19.1 

2.325 

.86 

120 

10.9 

1.5 

•  i 

.862 

2.95 

14-2 

21.7 

2.86 

.97 

125 

12,2 

2.0 

.  i 

.932 

3.51 

146 

23.3 

1.10 

129 

13.1 

2.5 

«  > 

.967 

4.70 

147 

24.6 

3.59 

1.29 

132 

14.1 

3.0 

.  ■ 

1.03 

4.60 

1  '■,;> 

1     *  / 

30.1 

3.92 

1.22 

134 

19.1 

1.5 

.3 

1.16 

1.4-7 

\55 

8.3 

3.31 

.62 

129 

1.2 

3.0 

1. 

2.10 

166 

10.5 

4.72 

.74 

140 

2.3 

1.5 

.5 

1.12 

1.08 

153 

4.0 

2.36 

.57 

125 

-1.3 

2.0 

.5 

1.21 

1.26 

156 

5,6 

3.40 

.55 

130 

-1.3 

2.5 

.5 

1.35 

1. 

161 

5.3 

3.99 

.60 

135 

-1.1 

2.0 

.7 

.932 

1.11 

147 

5.7 

2.415 

.57 

121 

-  .1 

2.0 

.9 

.72 

1.16 

122 

11.7 

.967 

.62 

103 

6.6 

Table  VT.   Data  obtained  in  experimental  determination  of  describing 
function  for  ramp  input.  See  figure  4.7  for  plot. 
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Another  interesting  fact  observed  from  the  plots  in  figure  4.7  is 

symmetry  a xmt  the  value  ~JL  =  0,5.  Thus  the  curve  for  £L  =  0.4 

TKv2  YKv2 

coincides  with  that  for  0,6,  the  curve  for  0,3  coincides  ■with  that  for 

0.7  etc, 

4.4  Verification  of  the  Proposed  Describing;  Function 

Although  the  experimental  results  presented  in  section  4.3  con- 
firmed the  existence  and  to  some  degree  the  nature  of  the  describing 
function,  another  experiment  was  designed  to  explicitly  verify  the 
theoretically  derived  describing  function.  In  addition  to  the  second 
and  third  order  control  systems  previously  simulated,  a,  fourth  order 

system  with  transfer  function  IL. ,  was  used.  In  con- 

S(S  +  1)(S  +  2)(S  +  4) 

ducting  these  runs  the  gain  of  each  linear  system  was  held  constant  so 

that  the  linear  curves  could  be  superimposed  on  the  plot  of  describing 

function  curves,  and  the  predicted  limit  cycles  were  obtained  from  the 

intersections.  The  value  of  the  ramp  input  was  varied  in  order  to  obtain 

different  values  of  the  ratio  ,  £,„.„  .  The  hysteresis  was  set  at  approxi- 

YK 

mately  one  volt  and  held  constant  throughout.  Nichols  Plots  of  the 
three  linear  systems  are  shown  in  figure  4.6  and  the  nondimensionalized 
describing  function  curves  are  shown  in  figure  4.2.  The  predicted  and 
observed  limit  cycles  for  each  system  are  compared  in  table  VII,  The 
Brush  recordings  are  displayed  in  figures  4,9  through  4,21. 

Since  some  of  the  predicted  limit  cycle  values  showed  dis- 
crepancies of  25$  or  more  from  those  observed,  it  is  felt  that  a  dis- 
cussion of  error  sources  is  warranted  here.  First  of  all  the  inherent 
errors  of  any  describing  function  analysis  are  present,  that  is  the 
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Fourth  Order  System 


B  . 

illy 

3 

X 

Recorded 
Freq, 

X 

Predicted 
Freq. 

;»Error 

A 

$Erro 
Freq. 

1 

.05 

.281 

3.16 

.209 

3.29 

.150 

3.9 

23.9 

2 

.10 

.562 

3.00 

.366 

2.18 

.44 

27.3 

20.2 

3 

.20 

1.125 

2.89 

.61 

2.12 

.75 

25.3 

23.0 

4 

.30 

1.687 

2.79 

.7*8 

2.35 

.87 

15.7 

8.2 

5 

.50 

2.812 

2.7^ 

.350 

2.56 

.93 

7.6 

9.4 

6 

.70 

3.937 

2.75 

.722 

2.35 

.37 

14.5 

21.0 

7 

•  GO 

4.400 

.572 

2.12 

.75 

22.0 

30.0 

Third  Order  ' 

System 

1 

.05 

.231 

1.74 

.379 

3.21 

.170 

84.0 

55.0 

2 

.10 

.562 

1.60 

.679 

1.77 

.56 

10,6 

17.5 

3 

.20 

1.125 

1.49 

1.17 

1.32 

1.22 

11.4 

4.3 

4 

.30 

1.637 

1.40 

1.48 

1.32 

1.60 

3.6 

8.1 

5 

.50 

2.812 

1.36 

1.6? 

1.40 

1.35 

2.9 

10.4 

6 

.70 

3.937 

1.45 

1.39 

1.32 

1.60 

9.0 

14.3 

7 

.80 

4.400 

1.53 

1.07 

1.32 

1.22 

13.7 

14.0 

Second  Order 

System 

1 

.05 

1.0 

1.21 

1.965 

2.52 

.75 

108 

62.0 

2 

.10 

2.0 

1.10 

3.31 

1.56 

2.4 

41.8 

27.5 

3 

.30 

1.5 

.62 

3.31 

.580 

3.37 

6.4 

1.8 

4 

.50 

1.5 

.57 

2.86 

.^90 

3.00 

14.0 

4.9 

5 

.70 

2.0 

.57 

2.41 

.504 

2.40 

11.5 

0.4 

o 

.90 

2.0 

.62 

.967 

1.25 

0.34 

101.5 

65.0 

Table  VII.   Recorded  and  predicted  limit  cycle  values  for  ramp  inputs. 
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describing  function  predicts  the  ratio  of  the  fundamental  sinusoidal 
components  of  the  input  and  output  signals.  In  many  of  the  runs  con- 
ducted (large  and  small  ratios  of  -~  )  the  input  to  the  non-linear 
element  was  decidely  non-sinuscicUl.  This  can  be  appreciated  by 
exaraining  the  Brush  recordings.  It  should  be  noted  that  the  value 
recorded  for  the  experimental  value  of  X  was  simply  ■§■  the  peak  to  peak 
amplitude  of  the  relay  input,  while  the  value  of  X  predicted  by  the 
describing  function  is  the  amplitude  of  the  fundamental  frequency 
coiroonent.  Appendix  II  shows  that  for  values  of  -JL.  near  0.1  or  0.9, 
and  for  values  of  ft  near  90°,  the  pull-in  angle  approaches  130°,  thus 
requiring  an  extremely  non-sinusoidal  input  signal  with  a  peak  to  peak 
a  rplitude  which  approaches  \   the  amplitude  of  the  fundamental  frequency 
co  rponent  as  illustrated  by  figure  h. 5. 

In  section  k,Z   and  appendix  II  the  limits  for  sinusoidal  opera- 
tion are  discussed  and  the  area  defining  these  limits  is  shown  on  the 
describing  function  curves  in  figure  4.2.  It  can  be  seen  from  the  table 
of  predicted  and  observed  values  that  the  errors  are  largest  outside 
the  area  of  sinusoidal  operation,  and  that  the  predicted  values  become 
increasingly  more  accurate  as  the  area  of  sinusoidal  operation  is  ap- 
proached. 

Another  consideration  in  evaluating  the  discrepancies  between  pre. 
dieted  and  observed  values  is  the  sensitivity  of  the  describing  function 
to  hysteresis.  The  analog  simulation  of  the  relay  employs  two  d.c. 

lifiers  in  the  analog  computer,  set  to  high  gain  values.  It  was  not 
uncommon  for  the  drift  in  these  amplifiers  to  cause  changes  in  the 
hysteresis  setting  of  several  tenths  of  a  volt.  In  order  to  gain  an 
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appreciation  for  nitude  of  errors  introduced  by  s  all  changes  in 

hysteresis,  a  representative  system  containing  fourth  order  linear  com- 
ponents was  examined.  The  results  of  this  examination  are  shoxm  in 
table  VIII  and  graphically  displayed  in  figure  4.8.  The  extreme  sensi- 
tivity of  the  amplitude  to  hysteresis  changes  is  evident.  It  is  also 

noteworthy  that  these  errors  are  essentially  invariant  with  the  ratio 

B 

- 

■2 

In  conclusion,  after  careful  analysis  of  the  data  recorded,  it  is 

apparent  that  the  general  validity  of  the  proposed  describing  function  is 

proven.  Specifically,  when  intersection  with  a  describing  function  curve 

occurs  within  the  area  defined  for  sinusoidal  operation  the  amplitude 

and  frequency  of  the  limit  cycle  are  predictable  within  10$  error,  '.hen 

the  intersection  occurs  outside  of  this  area  proportionately  greater 

errors  are  expected,  but  reasonable  and  useful  estimates  can  be  obtained. 

4.5  .■.'he  ?ia..rj  Describing  Function  for  Two-sided  Operation. 

If  the  sinusoidal  component  of  the  relay  input  is  greater  than 

the  critical  value  of  X.  defined  by  equation  4.12  or  4.19,  the  system 

will  enter  a  limit  cycle  in  which  both  sides  of  the  relay  are  actuated 

during  each  cycle.  In  this  case,  the  pull-in  and  drop-out  angles  must 

satisfy  equation  4.20,  which  is  obtained  from  the  simultaneous  solution 

of  equations  2.20  and  1-12. 

.22   <*9  +/92  -o(1  -^  =  2*JL 

'v2 

...  order  to  evaluate  the  describing  function  for  a  relay  with 

hysteresis  it  is  necessary  to  obtain  the  simultaneous  solution  of  equa- 
tion 4.22  and  equations  1-8  through  1-11.  This  solution  can  be  obtained 
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1. 

2 

,^Error 
W 

$3rror 

A 

^Error 

•ror 
X 

fMrror 

■  T 

&Error 

X 

Error 

r 

'Error 
X 

.3 

-30.8 

J±.6 

-16,3 

—  1  e  2 

+15.7 

+1.2 

+33.6 

+2,3 

A 

-30.2 

-2.2 

-16.3 

-1.1 

+15.9 

+1.1 

+35.1 

+3.3 

.5 

-31.3 

-16,1 

-2.2 

+15.7 

+1.1 

+35*3 

+2.2 

.6 

-30.2 

-2.2 

-16.3 

-1.1 

+15.9 

+1.1 

+35.1 

+3.3 

.7 

-30.3 

J^.6 

-1^.3 

-1.2 

+15.7 

+1.2 

+33.6 

+2.3 

Table  VIH.   Errors  in  predicted  amplitude  and  frequency  resulting  from 
small  relay  hysteresis  errors .  Zero  error  condition  is 
for  p-d  =  1.0.  Data  shown  is  for  a  system  with  fourth 
order  linear  components: 

2 

3(3  +  1)  (3  +  2)  (S  +  4)"~ 
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graphical:  o(  -ft     .  o(       -  fi .  for  constant  values 

of  X*  The  value  c  W  can  •    be  obtained  for  any  combination  of 

' •     and  X0  This  infoi  in  turn  be  used  to  calculate  de~ 

YK 
vp 

scribing  function  curves  which  will  be  valid  for  two-sided  operation, 
Unfortunately,  this  method  is  not  easily  adapted  to  a  non-dimension- 
alized  form.  Consequently  it  must  normally  be  calculated  individually 
for  each  relay  to  be  investigated. 

The  solution  is  slightly  less  complicated  for  a  relay  without 
hysteresis c  The  transfer  characteristics  for  such  a  relay  are  shown  in 
figure  3«1«  In  this  case,  of*  =  /3,  and  °^2  ~  ^2°  Therefore  the  system 
of  simultaneous  equations  which  must  be  solved  in  order  to  evaluate  W 
consists  of  only  equation  4,23  and  equations  1-9  and  1-11, 

The  equations  for  the  coefficients  of  the  fundamental  components 
of  the  fourier  series  expansion  of  0(t)  reduce  to 

4,24   a1  =  0 

4-0 25     -  ££  (cos  $  4   +  cos^«) 
IT  c 

The  solution  can  be  simplified  further  by  noting  that 

4,26       cos   $  os  fi  2  -  2  cosK^      -  fi  2)cos|(^1  -  $z) 


=      cosi(£«  +  /£?)cos       BTr 
d  2YK„ 


*.27      b1  =  iilcos   ( ^.    +  £    )cos  -JET 
^  2YK 


V 
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An  even  si  ion  can  be  obtained  for  a  system  containing 

an  ideal  relay.  In   this  case  there  is  no  dead  zone,  and  "^   -  -  & 2* 
Therefore  equation  4.2?  reduces  to  equation  4.28  and  the  describing 
function  can  be  expressed  by  equation  4.29. 


4.28   bf  =  i£L  cos 

77 


■TT 


It  should  be  noted  that  there  will  be  a  discontinuity  between  the 
describing  function  for  one-sided  operation  and  the  describing  function 
for  two-sided  operation  whenever  there  is  a  significant  hysteresis  effect. 
Therefore,  if  the  system  is  operating  on  only  one  side,  and  a  small 
change  in  input  or  gain  causes  the  error  signal  to  change  so  that 
X  >  Xc,  the  system  output  may  change  radically. 

It  is  possible  that  this  discontinuity  may  help  to  predict  some 

types  of  "jump  resonance"  found  in  relay  servos.  Suppose,  for  instance, 

that  the  system  is  subjected  to  a  sinusoidal  incut,  Q-    =   Asin  u>  t,  such 

that  ku  t  =  ,9¥Kv?«  ^  UQ   3-s  ^uch  less  than  the  natural  frequencies 

of  the  system,  the  servo  will  behave  as  if  the  input  signal  were  a 

slowly  changing  ramp  which  has  a  maximum  slope  of  Au>Q  at  u>Qt  =  0  and 

decreases  to  zero  at  \jJA>  -  ^   . 

2 

ow  suppose  the  system  characteristics  are  such  that  the  relay 
will  operate  on  only  one  side  when  <->0t  =  0e  Then  the  initial  operating 
point  will  be  determined  by  the  ramp  describing  function  for  one-sided 
operation.  The  operating  point  will  tend  to  travel  along  the  linear 
curve,  with  the  instantaneous  operating  point  being  determined  by  setting 
B  =  Ato»0cos  u>0te 
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If  the  linear  curve  crosses  the  critical  curve  of  D^  for  the 
relay  used,  the  operating  point  will  eventually  reach  a  point  where  one- 
sided operation  can  no  longer  occur,  and  a  two-sided  limit  cycle  will 
result.  During  the  next  portion  of  the  cycle,  the  servo  will  operate 
in  accordance  with  the  describing  function  for  two-sided  operation.  The 
servo  will  return  to  one-sided  operation  when  X  decreases  so  that 

Du  2E  ■ £——  (1  +  sin  o(U ) ,  where  both  X  and  ex.  ~  are  determined  by  the 

Pn  -  ^ 

tx;o-sided  describing  function. 

In  generals  the  magnitude  of  the  system  output  will  "jump"  to 
a  larger  value  when  the  relay  shifts  from  one-sided  operation  to  two- 
sided  operation,  and  will  "jump"  to  a  smaller  value  xjhenever  the  relay 
returns  to  one-sided  operation,  due  to  the  discontinuity  between  the  two 
describing  functions.,  There  may  also  be  a  change  in  the  frequency  of 
the  limit  cycle  when  the  shift  occurs „  Both  the  frequency  change  and 
magnitude  change  can  easily  be  predicted  once  the  describing  function 
curves  have  been  plotted 0 

^•.6  Extension  of  the  Ramp  Describing  Function  to  Other  Types  of  I'on- 
linearities., 

The  ramp  describing  function  can  be  extended  to  apply  to  other 
types  of  non-linearities.  In  general,  the  fourier  series  expansion  of 
the  output  wave  will  give  several  equations  of  the  following  types; 

1.  A  relationship  for  a0  similar  to  equation  1-3. 

2.  Relationships  for  a«  and  b*  similar  to  equations  1-5  and  1-7. 

3.  Relationships  for  all  angles  involved,  similar  to  equations 
I-S  through  1-1 
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Since  e  tes  for  any  non-linearity,  it  can  be 

combined  with  the  fourier  equivalent  of  a0  to  get  a  relation  between 
the  ramp  magnitude  and  the  angles  involved,  similar  to  equation  ^aLl-„ 
This  in  turn  can  be  combined  with  the  angle  relationships  to  solve  for 
the  dec«  component,  W.  Once  W  is  evaluated,  it  is  a  relatively  simple 
matter  to  evaluate  the  equations  for  a-j  and  b-j  and  substitute  these 
values  into  equation  3«3  to  obtain  the  describing  function* 

It  should  be  noted,  however,  that  the  transcendental  equations 
involved  will  become  more  complex  as  the  nature  of  the  non-linearity 
becomes  more  complex,  and  the  solution  of  these  equations  may  become  ex- 
tremely difficult,,  In  addition,  other  types  of  non-linearities  may  have 
different  modes  of  operation  corresponding  to  the  one-sided  and  two- 
sided  modes  of  operation  for  the  relay  servomechanism,  and  there  may  be 
discontinuities  between  the  describing  functions  for  these  different 
modes. 

Thus  a  complete  analysis  will  require  the  evaluation  of  a  separate 
group  of  describing  function  curves  for  each  mode  of  operation,  along 
with  a  set  of  limit  curves  to  define  the  points  at  which  the  system  will 
shift  from  one  mode  of  ooe ration  to  another. 
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.  ft  V 
CONS  FOR  FUTURE  INSTIGATIONS 

During  the  course  of  this  investigation,  several  interesting 
problems  were  noted  which  appeared  to  be  worthy  of  further  study.  These 
suggested  subjects  are: 

1 .  A  comprehensive  study  of  the  effect  of  different  types  of 
dissymmetry  on  systems  containing  different  types  of  non- 
linearities  . 

2.  A  more  exhaustive  analysis  of  the  ramp  describing  function 
for  a  relay  servo-mechanism  in  which  the  relay  operates  on 
both  sides.  This  analysis  should  include  an  investigation 
of  the  discontinuity  between  the  describing  function  for  one- 
sided operation  and  the  describing  function  for  two-sided 
operation  to  determine  whether  or  not  these  describing 
functions  can  be  used  to  predict  certain  types  of  "jump 
resonance"  effects,  as  suggested  in  section  4.5. 

3.  A  thorough  study  of  the  extension  of  the  ramp  describing 
function  to  non-linear  devices  other  than  the  relay  would  also 
be  valuable  to  determine  the  effect  of  ramp  inputs  on  systems 
containing  such  elements.  An  outline  for  the  procedure  for 
such  an  investigation  is  given  in  section  4.6. 
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.  APPEKDIX  I 

DERIVATION  OF  THE  SSi^EfULIZED  DESCRIBING  FUNCTION 

FOR  A  DI33Y2  METRIC  RELAY 


A.     Input  vs„   output  characteristics. 


wt  *  7T  -  /3, 


cjt'ir*  *(. 


UJt*2fT-&2 


Figure  1=1 
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0(t)  =  X., 
0(t)  =  0 
0(t)  =  Y2 
.:.)  =  o 
Fourier  series  expansion  of  output  irave 


0  <  ut    -    °^| 

o<1  <  tot      <     77     -  #j 
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0(t)  =  _2     +  a.*   cos  cot  +  b-j    sin  u)t 
0(t)  =  %■  4  Csin(u>t  +  0) 

c=  yV  •  bi2 

0  =  tan"1     fl 
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a.  =  J- 
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0(t)  du>t  =  '. 
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TT    +  crt 
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1-3     a0  =  l         ;tr.  /9,  -  «,)  -  V^-^2-^2) 


'2TT 
a1         -  I      0(t)   cosuit  dwt 


'77-/8 


i^ cos  tot  d  wt  +  1 
TT 


2TT„0, 


-^cosut  dwt 


77 


Y> 
^  -  sino^)  +    ~(s2iifi2  -  sino<9) 

,   =  .-  /     O(t)   sin  uj  t    du)t 

1    A77  ^2 

sin  uJt  d  *^t  +  —  /         r0sini*>t  d«Jt 

Vrr  +o<9 

Y  Y 

1-7     b1   ■:  ~^(cos  o^        cos  /3   )    f  =r(cos  o(2  +  cos  $J 

For  a  perfect  sine  input,   the  angles  can  be  defined  as  follows: 

1-8       sin  o^   =  ~^~  =  Rpi 

n         d1    "   W 
1.9       sin^  =,-!_    =Rdi 

P9    -    VJ 

1-10     sin  o(?  =   ^-J-    =  ^p2 

/?        do  -  W 
1-11   sin  P2  =  -2 = 

C.     The  describing  function. 

The  generalized  describing  function  is: 


3.3       Gd=l   J**  +  b,2  /,n-a     g 


b1 

If  both  sices  of  the  relay  are  actuated  on  each  cycle ,  equations 
1-4  through  1-7  can  be  used  in  the  evaluation  of  equation  3»8  without  re. 
strictions.  If  the  relay  is  only  partially  actuated,  the  same  equations 
still  be  used,  pro      certain  angles  are  defined  as  follows: 

If  -o  >    1,      o<  =  ?  n  =  1,2 
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If  -V   <  -1,  then  °*.  =  =?  n  =  1,2 

—  -n   ^ 

If  Hrf   >  1 ,  then  £  =  ~  n  =  1 , 2 

^n  n   ^ 

If  R^  >  -1,  then  ^n==f  n  =  1,2 

D,  Special  Cases  for  "^  =  1^ 

If  I*  =  lo  =  I,  equations  1-3,  1-5 »  and  1-7  reduce  to 

1-12   aQ  =  fcfiz  +  ^2  -  *,  -  *,) 

1-1 3   a -j  =  ~(sin  /5  ^  -  sin  o<   +  sin/?  -  sin  <^  ) 

1-14   b^  =  ^(cos  j3   -j  +   cos  o(  .  +  cos  /6  2  +  cos  o<  ) 

If  ir-j  =  Y2  =  Y  and  the  relay  has  no  hysteresis,  then  °( *   =  8 * 
and  <^2  =  ^2»  anc*  equations  1-12,  1-13  and  1-1' l  reduce  to 

1-15   ao=f(^2-^l) 
3.6    a1  =  0 


1-16   b-j  =  ~(cos  /S,  +  cos  ySJ 


t^>  v*   v   —  V 


.j  =  Y«  =  ^  and  the  relay  operates  only  on  the  positive  side, 

then  P  o  =  o(  ~  -    —  e   The  subscripts  can  be  dromed  for  convenience  in 
^     -    2 

writing,  so  that  equations  1-12,  1-1 3,  and  1-14  reduce  to 

4.1  a0=£(7r-£-O 

4.2  a^  =  l(sin  /S  -  sin  c<  ) 

1   /^ 
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TT 


(cos  £  +  cos  o(  ) 
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rSOIDAL  OPERATION  DURING 

:  . :  3iDS 

relay  operates  on  only  one  side,  sinusoidal  operation- 
is  possible  only  if  c<  —  90°. 

*A      ^2-=  (77-  /3  .  <<) 

and 


Therefore, 


*.7   £  -  ^-^ 


U-1  4=2"-  ^JL    +  0 

-v2 


<*  =  II  -  BIT 

Cv2 


2   _, 


For  bhe  case  where  /6  <  -  —  ,  equation  Ii-1  gives 

2 

n-3   -t  >  rf  -  -£5L 

^2 


srefore.    if        >:    ,     TLT  ,    sinusoidal  operation  is  possible  only  if 

^2 

^.20   -£  ^  77-  -231 

v2 

Tor  the  case  where  o^  >  H  ,  equation  II-2  reduces  to 

:v2 
refore  if  ]   ^  ,5IKV  ,  sinusoidal  operation  is  possible  only  if 


•V2! 


4.21   -jfi  <  «22L 


3ir 


-v2 
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